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43 \ Radiatively generated, LO quark (u, d, s, c, b) and gluon densities in the real, unpolar- 

ized photon, improved in respect to our previous paper [1], are presented. We perform 
three global fits to the F^ data, using the LO DGLAP evolution equation. We improve 
the treatment of the strong coupling running and used lower values of Aqcd, as we have 
found that the too high values adopted in the previous work caused the high x 2 of the fits. 
In addition to the modified FFNSc* jkl model, referred to as FFNScj^-1 we analyse a 
FFNScj^-2 model in which we take into account the resolved-photon heavy-quark contri- 
bution. New CJK model with an improved high-x behaviour of the F^(x, Q 2 ) is proposed. 
Finally, in the case of the CJK model we abandon the valence sum rule imposed on the 
VMD input densities. New fits give x 2 P er degree of freedom about 0.25 better than the 
old results. All features of the CJKL model, such as the realistic heavy-quark distribu- 
tions, good description of the LEP data on the Q 2 dependence of the F^ and on F^ c are 
preserved. Moreover we present results of an analysis of the uncertainties of the CJK par- 
ton distributions due to the experimental errors. It is based on the Hessian method used 
for the proton and very recently applied for the photon by one of us. Parton and struc- 
ture function parametrizations of the best fits in both FFNSc*jj<- and CJK approaches are 
made accessible. For the CJK model we provide also sets of test parametrizations which 
allow for calculation of uncertainties of any physical value depending on the real photon 
parton densities. 



1 



1 Introduction 



In this paper we continue our recent analysis of the LO unpolarized real photon parton distri- 
butions, [1], improving and broadening our research. The main topic of our previous work was 
the description, within the DGLAP evolution framework, of the heavy, charm- and bottom-, 
quark contributions to the photon structure-function, F^(x, Q 2 ). We analysed and compared 
two approaches In the first analysed model, referred to as FFNScjkl, we adopted a widely 
used massive quark approach in which heavy quark, h, contributes to the photon structure 
only through the so-called Bethe-Heitler, 7*7 — > hh process. In such models the heavy-quark 
masses are kept to their physical values. In the second, CJKL model, we used the ACOT(x) [4] 
scheme, where heavy-quark densities appear. It was the very first application of this scheme to 
the photon structure. We performed two global fits to the set of updated F^(x, Q 2 ) data col- 
lected in various e + e _ experiments. We based both models on the idea of radiatively generated 
parton distributions introduced by the GRV group (see [5] for the photonic case). 

In this work the main assumptions of the previous analysis are left unchanged, although 
some details are improved, and we analyse an additional model. First of all, we improve the 
description of the running of the strong coupling constant, a s , and use Aq CD v& hie substantially 
smaller than the one used previously. The Aq C d values applied in the former analysis were 
obtained with the assumption that the LO and NLO Aq C d values for four active flavours are 
equal, as in the GRV analysis. They appeared to be too high and were one of the causes for 
the high \ 2 of the fits. 

The old FFNScjxl model is now being denoted as FFNScxkI- We compare it with the more 
realistic model FFNScj/^, in which we include the so-called "resolved-photon" contribution 
to F^x^Q 2 ) given by the 7*G 7 — > hh process [6]. As we stated already in [1], the CJKL type 
model needs correction which could resolve the problem of the high-x behaviour of F^{x, Q 2 ), 
predicted by this model. Therefore, we analyse a new CJK model in which this problem 
is avoided. The CJK model is further slightly improved by changing the lower limit of the 
integration, in the heavy-quark subtraction terms, from the square of the mass of the heavy 
quark, m^, to the starting scale of the DGLAP evolution, Q%. This also leads to better fits. 
Finally in the case of the CJK model we abandon the valence-number sum rule imposed by 
hand on the VMD input densities but keep the corresponding energy- momentum constraint. 

Finally, in all new fits we use limited set of data, excluding the TPC27 data which are 
considered to be inconsistent with other experimental results. This is an additional reason, 
beyond improvements mentioned above, why our new predictions for model FFNScjif 1 are not 
identical with results obtained in [1] from full set of data. 

We also proceed in the direction not addressed before. We analyse uncertainties of the 
parton distributions due to the experimental errors of F^(x, Q 2 ) data. This part of the work 
has been motivated by the recent analysis performed for the proton structure by the CTEQ 
Collaboration, [7-9] and the MRST group, [10]. We use the Hessian method, formulated in 
recent papers, to obtain sets of test parton densities allowing along with the parton distributions 
of the best fit to calculate the best estimate and uncertainty of any observable depending on 
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the photon structure. Full discussion has been given in [3], see also [2]. 

This paper is divided into four parts. Section 2 recalls the previous FFNScxkx and CJKL 
models of the real photon structure. In section 3 we explain in detail all the changes currently 
introduced to the models. Next, in section 4, we present results of the new fits to the experi- 
mental F^(x, Q 2 ) data along with the calculated uncertainties of the CJK parton distributions. 
The parton distributions which are a result of our analysis have been parametrized on the 
grid. In section 5 we give a short summary and information where to find the corresponding 
FORTAN programs. 



2 FFNS^j^l and CJKL models - short recollection 

The two approaches leading to the FFNSc jkl and CJKL models, and considered also in this 
analysis, have been described in detail in our previous paper [1]. The difference between them 
lies in the approach to the calculation of the heavy, charm- and beauty-, quark contributions to 
the photon structure function F^x.Q 2 ). First, FFNScja-l model bases on a widely adopted 
Fixed Flavour-Number Scheme in which there are no heavy quarks (denoted below by h) as 
partons in the photon. Their contributions to F^(x, Q 2 ) are given by the 'direct' (Bethe-Heitler) 
7*7 — > hh process. In addition one can also include the so-called 'resolved'-photon contribution: 
7*G 7 — > hh. We denote these terms as F^ h {x, Q 2 )\di r and F^ h (x, Q 2 )\ res , respectively. In this 
paper we consider two FFNS models: in the first one, FFNScjatI, we neglect the resolved- 
photon contribution, while in the second one, FFNScja-2, both mentioned contributions to 
Fi](x,Q 2 ) are included. The photon structure function is then computed as 

3 

F;\x,Q 2 ) =$> ei 2 (g7 + q?)(x,Q 2 ) + [F2, h (x,Q 2 )\ dir + F2 ih (x,Q 2 )\ res ], (1) 

i=l h(=c,b) 

with q^(x,Q 2 ) (qJ(x,Q 2 )) being the light u,d,s quark (anti-quark) densities, governed by the 
LO Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations, [11]. 

The CJKL model adopts the new ACOT(x) scheme, [4], which is a recent realization of the 
General Variable-Flavour Number Scheme (GVFNS). In this scheme one combines the Zero- 
mass Variable-Flavour Number Scheme (ZVFNS), where the heavy quarks are considered as 
massless partons of the photon, with the FFNS just discussed above. In this model, in addition 
to the terms shown in Eq. (0), one must include the contributions due to the heavy-quark 
densities which now appear also in the DGLAP evolution equations. A double counting of the 
heavy-quark contributions to F%(x, Q 2 ) must be corrected with the introduction of subtraction 
terms for both, the direct- and resolved-photon, contributions. Further, following the ACOT(x) 
scheme, we introduce the Xh — x(l+4m^/ Q 2 ) variables giving the proper vanishing of the heavy- 
quark densities at the kinematic thresholds for their production in DIS: W 2 = Q 2 (l — x)/x > 
4m^, where W is the 7*7 centre of mass energy. Adequate substitution of x with \h i n Qh and 
the subtraction terms forces their correct threshold behaviour, as Xh — > 1 when W — > 2rrih- 
This is achieved for all terms except for the direct subtraction term FiJ^x, Q 2 )\dir,sub for which 
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there is a need of an additional condition, F^ h (x, Q 2 )\dir,sub = for \h > 1- The full formula 
for the photon structure function in the CJKL model is 

F?(x,Q 2 ) = xj^e^ql + qlMx^ + x £ e\W h + <f h )( X H, Q 2 ) 

i=l h{=c,b) 

+ [F^ h (x,Q 2 )\ dvr + Fl h (x,Q 2 )\ res ] (2) 

h(=c,b) 

~ ^ 4 [F^hldir^ubiXhj Q ) + F^ h \res,sub{Xh, Q )] • 
h{=c,b) 

with positivity constraint for each heavy-quark contribution , F^^x^Q 2 ) > 0. Explicit expres- 
sions for the terms appearing in Eqs. and (J2J) can be found in [1]. 

We use the DGLAP equations summing the QCD corrections in form of leading logarithms 
of Q 2 . Their solution for quark densities can be divided into the so-called point-like (pi) part, 
equal to a special solution of the full inhomogeneous equation and the hadron-like (had) part, 
arising as a general solution of the homogeneous equation. Their sum gives the partonic density 
in the photon. It is most useful to write the result in the Mellin-moments space with the nth 
moment defined as ^ 

f n (Q 2 )= [ x^ffaQ^dx, (3) 
Jo 

where f(x,Q 2 ) can be the parton (quark and gluon) densities, g 7 (x, Q 2 ), or the splitting func- 
tions in the DGLAP evolution equations P(x, Q 2 ) and k(x, Q 2 ) defined in [1]. Then one obtains 



where 



(5) 

&Q 2 ) = L~ 2pn '^{Ql). 

Here L = ""^j , where is the scale at which the evolution starts (we call it the input scale). 

For all models we choose to start the DGLAP evolution at small value of the Q 2 scale, 
Qq = 0.25 GeV 2 , following GRV [5] hence our parton densities are radiatively generated. As 
it is seen above the point-like contributions are calculable without further assumptions, while 
the hadronic parts need input distributions. For this purpose we utilize the Vector Meson 
Dominance (VMD) model [12], where 

fUx,Qt) = T,^f v ( x >Q°)> ( 6 ) 

v Jv 

with the sum running over all light vector mesons (V) into which the photon can fluctuate. 
The parameters f v can be extracted from the experimental data on the T(V — *> e + e~) width. 
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In this analysis we use explicitely the p°-meson densities while the contributions from other 
mesons are accounted for via a parameter k, that is left as a free parameter. Thus, we take the 
parton densities in the photon equal to 

fL d ^Ql) = ^f p ^QD- (7) 
fp 

We take the input densities of the p° meson at Ql = 0.25 GeV 2 in the form of valence-like 
distributions both for the (light) quark {v p ) and gluon (G p ) densities. All sea-quark distributions 
(denoted as ( p ), including s-quarks, are neglected at the input scale. 

The valence-quark and gluon densities satisfy the energy- momentum sum rule for p°: 

x(2v p (x,Ql) + G p (x,Ql))dx = 1, (8) 

and the sum rule related to the number of valence quarks, n v 

n v = [ 2v p (x,Q 2 )=2. (9) 
Jo 

Both of them we imposed as constraints on the parameters of the models in our previous 
analysis [1]. 

The input quark and gluon densities are taken in the form (with a > 0) 

xv p (x,Ql) = N v x a (l-xf, 

xG p (x,Ql) = N g xv p (x,Q 2 ) = N g x a (l-x)P, (10) 
x( p (x,Ql) = 0, 

where N g = N g N v . The imposed constraint given by Eq. (jHJ) allows to express the normalization 
factor N g as a function of a,/3 and N v . Moreover, when constraint is imposed the N v 
parameter can be further expressed in terms of a and (3. In the former case there are three and 
in the letter case four free parameters left as a subject to the global fit to the F^x, Q 2 ) data. 



3 New analysis 

We have performed new fits with a slightly changed data set as compared to the previous work. 
Moreover, in our new analysis we improved the treatment of the running of a s , by differentiating 
the number of active quarks in the running of a s and in the evolution equations, and by using 
lower values of Aqcd- We first describe new aspects of our analysis which are common to all 
considered models. Aspects relevant only for the CJK model are discussed next. 
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3.1 Data Set 



New fits were performed using all existing F^(x,Q 2 ) data, [14-24], apart from the old TPC27, 
[25]. In our former global analysis [1] we used 208 F^(x, Q 2 ) experimental points. Now we 
decided to exclude the TPC27 data from the set because it has been pointed out (see for 
instance [26]) that these data are not in agreement with other measurements. A detailed study 
of the influence of various F% data sets on the fits is performed by one of us in [3]. After the 
exclusion of the TPC27 data we are left with 182 F](x, Q 2 ) experimental points. We include 
all these data in the x 2 fit without any weights. A list of all experimental points used can be 
found on the web-page [27]. 

The exclusion of the TPC27 points affects the xVdof of the fits but has very small influence 
on the shape of the resulting parton distributions. 



3.2 a s running and values of A^) 

The running of the strong coupling constant at lowest order is given by the well-known formula: 

a '^ {Q2) = wk^r with a = (») 

where N q is the number of quarks entering in the a s evaluation. 1 This number increases by one 
unit whenever Q 2 reaches a heavy-quark threshold, i.e. when Q 2 = m\, where the condition 
ai Nq \m1) = o^ q+1 \mfy is imposed in order to ensure the continuity of the strong coupling 
constant. 

In our previous analysis Nj was identified with the number of active quarks in the photon: 
Nf = 3 and 5 in the FFNScjifL and CJKL models, respectively. Since we now distinguish 
between both numbers of quarks we have to use slightly more complicated formulea for the 
evolution of the parton densities, as now the above equations depend also on N q through 
their dependence on a s (Q 2 ) and Because of the implicit introduction of the heavy-quark 
thresholds into the a s running we must proceed in three steps to perform the DGLAP evolution. 
In the first step, describing the evolution from the input scale Q to the charm-quark mass m c , 
the hadronic input q^ ad {x, Qq) is taken from the VMD model. In the second step we evolve 
the parton distributions from m c to the beauty-quark mass, mi, and in the third one we start 
at mj,. In the second and third steps a new hadronic input at Q 2 = m\ is given by the sum 
of the already evolved hadronic and point-like contributions and the point-like distribution at 
this scale becomes zero again. 

In the previous work we assumed (following the GRV group approach [5]) that the LO and 
NLO Aqcd values for four active flavours are equal. We adopted Aq^ d = 280 MeV, value 
given in the Particle Data Group (PDG) report [13]. We now abandon this assumption and 

1 Notice that we distinguish now between the number of active quarks in the photon, denoted by Nf, and 
the number of quarks contributing to the running of a s , N q . 
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take Aq CD = 115 MeV, which is obtained from the world average value a s {Mz) = 0.117, with 
Mz = 91.188 GeV, using the LO expression for a s evolution, Eq. (JTTJ). As a consistency check 
we performed fits keeping JS^q CD as a free parameter and obtained results close to 115 MeV. 
Since it is not our aim in this paper to extract a value of a s from a fit to F% data, we prefer to 
fix A.q CD = 115 MeV rather than add a new free parameter in our fits. Imposing the continuity 
condition for the strong coupling constant and m c = 1.3 GeV and m& = 4.3 GeV, we obtain 
A { nl n = 138 MeV and A^L = 84 MeV. 



3.3 VMD input 

Like in our previous work the input evolution scale has been chosen to be small, Qq = 0.25 
GeV 2 for both types of models and we apply the same form of the VMD model input, given by 
Eqs. flmj) and Eq. ©. 

Finally, we try to relax the constraint on number of the valence quarks, n v , in the p meson. 
This leads to 4-parameter fits. That, as will be quantitatively shown in next sections, is possible 
only in the case of the CJK model. Therefore in each of the new FFNScxr- models we have 3 
free parameters. 



3.4 Modified subtraction terms in the CJK model 

In [1] we derived the subtraction term for a direct contribution, F% h \dir,sub, from the integration 
of a part of the DGLAP evolution equations, namely: 

where k(x) is the lowest order photon-quark splitting function (see Eq. (7) in [1]). The question 
here is: What should the limits of the integration be? The upper limit is obviously the Q 2 
scale at which the subtraction term is calculated. For the lower limit we took previously the 
standard for the Bethe-Heitler process scale: Qf ow = ml. However, the threshold condition is 
W 2 < 4m 2 h . This means that even for Q 2 < m 2 the heavy-quark contributions do not vanish as 
long as the condition x < Q 2 / (Q 2 + 4m^) is fulfilled. Hence, in this paper we take Qf ow = Ql 
and the direct subtraction term is given by: 

F2 h \dir,sub(x, Q 2 ) = x\n% I ?>e A h -(x 2 + (1 - x) 2 ). (13) 

Qo 71 

The same discussion applies to the subtraction term for the resolved-photon contribution. So 
we now use 

i^U,^ = x\n%e 2 h ^l f ^P qG (-)G^y,Q 2 ). (14) 

vo n j x y y 

instead of Eq. (19) in [1]. We found that the quality of the fit improves with this choice of the 

QloW 
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As we noticed in Section 2 the x — > Xh substitution leads to the proper threshold behaviour 
of all the heavy-quark contributions to the F^(x, Q 2 ), except for the subtraction term for direct 
contribution. It is already seen in Eq. (JT3Jl that this term does not vanish for Xh ~~ * 1 an d 
therefore by subtracting it the resulting heavy-quark contribution to may become negative 
in some regions of the x and Q 2 plane. An extra constraint to avoid this unphysical situation 
is, thus, needed. In Ref. [1] we imposed the simple condition (positivity constraint) that the 
heavy-quark contribution to F^{x, Q 2 ) has to be positive. Unfortunately, this constraint was 
not strong enough and for some small windows at small and large x still the unphysical situation 
^2hi x i Q 2 ) < ^2h( x y Q 2 )\dir + F2 h (x, Q 2 )\ re s was found [28]. Therefore, in this analysis we apply 
a positivity condition in the following form: 

Fl h {x,Q 2 ) > FZ h (x,Q 2 )\ dtr + FZ h (x,Q 2 )\ res . (15) 

4 Results of the new F 2 7 global fits 

In this analysis we determine the parameters of the models, related to the initial quark and 
gluon densities at the scale Ql = 0.25 GeV 2 , by means of the global fits to the experimental 
data on F^x^Q 2 ). We use 182 F^x^Q 2 ) experimental points, [14-24], with equal weights. 
Still, it has been shown in [3] that when we remove the CELLO [14] and DELPHI 2001 [20] 
data sets the quality of the fit improves substantially but the parton distributions lie well within 
the CJK uncertainties. However, we think that there is no strong argument to discard these 
data sets. So, we have chosen to keep them with the same weight as the other sets. Fits 
are based on the least-squares principle (minimum of x 2 ) an d were done using Minuit [29]. 
Systematic and statistical errors on data points were added in quadrature. 

In the CJK model we have four free parameters: a,(3,N v ,K, Eqs. (fTUJl and Eq. (JJJ). On 
the other hand, the two FFNS models differ only in the inclusion or not of the resolved- 
photon contribution to F^(x, Q 2 )(only the FFNScj^ model takes it into account through the 
7*G 7 — > hh process). For both FFNS models we impose the number of valence quarks constraint 
(J9~|) that allows to express N v in terms of a and (3 reducing the number of free parameters to 
three. 

The parameters of our new fits are presented in Table ^ The second and third columns 
show the quality of the fits, i.e. the total x 2 f° r 182 points and the x 2 P er degree of freedom. 
The fitted values for parameters a, f3, k and N v are presented in the middle of the table with 
the errors obtained from Minos with the standard requirement of Ax 2 = 1. In the last column 
the value obtained from the constraint (jSJ) for N g from other parameters is given. 

In the case of the FFNScjif models the test fits with the abandoned constraint Q gives 
n v ~ 0.5 and ~ 1.4 in the FFNScj/rl and FFNScjif2 models, respectively. We think that it is 
too far away from the expected value n v = 2 2 . This is the reason for keeping the constraint 
© for both FFNScxr- models. In this case the N v parameter is calculated from the constraint 
© and therefore we do not state its error. 

2 Note, that in the CJK model we obtain n v = 2.0 ± 0.1 
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model 


X 2 (182 pts) 


X 2 /dof 


K 


a 




N v 


N 9 


FFNScj^I 


314.0 


1.754 


9 o«7+0.063 
Z.ZU 1 _ .072 




f) 7Q9+0- 1 89 
u - ' yz -0.149 


0.358 


5.02 


FFNS C j*2 


279.8 


1.563 


o 1 i n +0.084 
Z.l±U_g 090 


n Qin+o.054 
u.oiu_ 051 


fl Q9Q+0.265 
U.OZO_ .223 


0.415 


4.51 


CJK 


273.7 


1.537 


I Q34+0.131 
1 - 3 °^0.124 


u.zyy_ 069 


0.898jg" 


D 4D4 +0 - 116 


4.93 



Table 1: Results of the fits for the three models considered in the analysis. The quoted errors 
are obtained from MINOS with the standard requirement of A% 2 = 1. 



The x 2 P er degree of freedom presented in Table ^ are still rather high, however there is 
an improvement as compared with previous fits. The old x 2 /dof for the same set of 182 data 
points read 1.99 in the FFNScjxl and 1.80 in the CJKL model. We see that the corresponding 
new xVdof values are about 0.25 lower. This is mostly due to the adoption of much lower 
Aqcd values as well as the modification of the subtraction contributions in the CJK model. 
The rejection of the TPC27 data also reduces the value of xVdof- 

We observe that the xVdof for the FFNScjif2 and CJK models are very similar and 
lower than the one for the FFNScj/^l model. It is obvious that the inclusion of the resolved 
7*G 7 — > hh contribution to F^(x, Q 2 )improves the agreement between the model and the data 
(see for instance [30]). 

We see that the quality of the present data does not allow for a clear discrimination between 
the different ways of dealing with the heavy quarks as the \ 2 and all fitted parameters are very 
similar. The k values are close to 2 which is in agreement with the GRV LO [5] prediction. 
The a parameter varies from about 0.25 to 0.31 (a — 1 ~ —0.75, —0.7) which seems to be 
conformable with the Regge theory predicting that for a valence-quark density a — 1 ~ —0.5. 
The (3 values ranging from about 0.8 to 0.9 are again in good agreement with the GRV LO [31] 
finding (0.85). They are not as small as in the case of our former FFNScxr-z, model but are far 
from 2, a standard prediction from the quark-counting rule [32]. 

4.1 Comparison of the CJK and FFNScjkI, 2 fits with the F 2 7 data 

In this paper we are going to present plots only for the three models that we analyse, without 
any comparison with other parametrizations. However, we will present plots of tthe same 
type as in Ref. [1] in order to facilitate the comparison. Moreover we will describe differences 
between of our new results and the previous ones and with the GRS LO [30] and SaSID [33] 
parametrizations. If experimental points for a few values of Q 2 are displayed in a panel, 
the average of the smallest and biggest one was taken in the computation of the theoretical 
prediction. 

Figures HHH show a comparison of the CJK and FFNScjif fits to the F^x^Q 2 ) with the 
experimental data as a function of x, for different values of Q 2 . The FFNScjifl fit predictions 
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are very similar to the GRS LO parametrization results in the whole range of x while the 
FFNScxr-2 and the CJK model predict a much steeper behaviour of the F^x^Q 2 ) at small x 
with respect to the FFNS^j^l fit (and GRS LO) and SaSID parametrizations. On the other 
hand these curves are less steep than the old FFNSc jkl and CJKL ones. The behaviour of 
the three fits in the region x > 0.1, as shown in Figs. EJandlH is very similar. 

Apart from this direct comparison with the photon structure function F% data, we perform 
another comparison, this time with LEP data that were not used directly in our analysis. 
Figures El and El present our predictions for the F^x.Q 2 ), averaged over various x regions, 
compared with the recent OPAL data [24]. Like in our previous analysis we see that all FFNS 
type predictions (including GRS LO and SaSID parametrizations) are similar and fairly well 
describe the experimental data. Moreover, again the CJK model, alike the CJKL model, slightly 
differs from other fits. However, this difference is much smaller now and gives better agreement 
with the data. 

We observe that for the case of the medium-x range, 0.1 < x < 0.6, there are small differ- 
ences between the CJK and both FFNS models at very low Q 2 , where there are no experimental 
data, and at large Q 2 where a slightly better agreement between the CJK model prediction 
and the OPAL data is found. Comparing with the plot in Fig. El we see that larger values of 
the F^(x, Q 2 ) obtained for the CJK model are originated at lower values of x in the considered 
range, as one could expect from Figs. HHU 

4.2 Parton densities 

It is very instructive to discuss the parton densities obtained in the CJK and FFNScxr- fits. 
We choose to present results for medium- and high-x at Q 2 = 10 GeV 2 , see Fig. In Figs. 
IHlfTTil we show the up-, and charm-quark and gluon densities for various Q 2 values. First we 
observe that parton distributions obtained in various models are very similar, except that of 
course there are no heavy-quark distributions in FFNS-type approaches. In the case of the 
CJK models, due to the introduction of the \h variable, the c 7 (x,Q 2 ) and fe 7 (x, Q 2 ) densities 
vanish not at x = 1, as in the case of the GRV LO [5] and SaSID parametrizations, but as 
they should at the kinematic threshold. Moreover, the CJK heavy-quark densities at low x are 
larger than the corresponding densities obtained in other parametrizations. This is a feature 
that can be observed in a wide range of Q 2 values in Fig. El 

We notice that our new parton densities have all similar shapes to the corresponding old 
CJKL distributions. Though, in the medium- and high-x regions they have slightly higher 
values. On contrary, at very low x values CJKL densities are greater than the new ones, 
opposite case occurs in the 10 -3 < x < 10~ 2 region. In the case of the gluon density we find 
that all new lines are much steeper at high-x than the predictions of our previous models and 
the GRV LO and SaSID parametrizations. 
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4.3 Comparison with F^ c 

We compare the individual contributions included in the C JK model relevant for our predictions 
of the F^ c . Results from the CJK fit are presented in Fig. HUfor Q 2 = 5, 20, 100 and 1000 GeV 2 . 
Almost all the contributing terms vanish in the W — > 2m c threshold in a natural way due to 
the introduction of the Xh variable. The only exception is the one subtraction term, namely 
F2, c \dir,sub which dominates near the highest kinematically allowed x value and vanishes only 
due to the extra condition, F2 )C \dir,sub = for Xh > 1, we imposed. The direct (Bethe- Heftier) 
term is important in the medium- x range. Its shape resembles the valence-type distribution. 
The charm-quark density contribution, i.e. the term 2xe 2 c l {x, Q 2 ), is important in the whole 
kinematically available x range and dominates the for small x. In this region also both 
resolved-photon contributions increase with decreasing x, but to great extent they cancel each 
other. 

Finally we see that imposing the improved positivity condition F% h > F 2: h(x, Q 2 ) 7 (x, Q 2 ) \di r + 
F^ h (x,Q 2 )\ res on the heavy-quark contributions to the F^(x,Q 2 ) results in correct threshold 
behaviour of the total i^c function. Unlike in the case of the CJKL fit, the Fl h and its 
contributions vanish at the same high x value. 

A good test of the charm-quark contributions is provided by the OPAL measurement of the 
F^ci obtained from the inclusive production of D* ± mesons in photon-photon collisions [34]. 
The averaged F^ has been determined in the two x bins. These data points are compared to 
the predictions of the CJK and FFNScxr- models and GRS LO and SaSID parametrizations 
in Fig. D21 

Our first observation is the following: our models containing the resolved-photon contribu- 
tion, F2 h (x, Q 2 )\res (FFNScxr:2 and the CJK model) agree better with the low-x experimental 
point than other predictions. The GRS LO and SaSID parametrizations also include the 
resolved-photon term but in their case the gluon density increased less steep than our models 
predict, as was already mentioned. Their lines lie below the results of our new fits but 
higher than the FFNScxr- 1 curve, given solely by the direct Bethe- Heitler contribution. 

The CJK model overshoots the experimental point at high x while other predictions agree 
with it within its uncertainty bounds. Taking into account both data points we conclude that 
the best agreement with the experimental results is provided by the FFNScxr:2 model. 

4.4 Gluon densities at HERA 

We also checked that the gluon densities of CJK and FFNScxk models agree with the HI 
measurement of the G 1 distribution performed at Q 2 = 74 GeV 2 [35]. As can be seen in 
Fig. Uni all models predict gluon densities that lie above the one provided by the GRV LO 
parametrization, which gave so far best agreement with the HI data. Further comparison of 
our gluon densities to the HI data cannot be performed in a fully consistent way, since the 
GRV LO proton and photon parametrization were used in the experiment in order to extract 
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such gluon density. 



4.5 The uncertainties of the CJK parton distributions 

Let us now present the main results of an analysis of the CJK parton distribution uncertainties 
described in detail in [3]. 

During the last two years numerous analysis of the uncertainties of the proton parton den- 
sities resulting from the experimental data errors appeared. The CTEQ Collaboration in a 
series of publications, [7-9], developed and applied a new method of their treatment signifi- 
cantly improving the traditional approach to this matter. Later the same approach has been 
applied by the MRST group in [10]. The method, denoted as Hessian method, bases on the 
Hessian formalism. We applied it for the very first time to the case of the photon parton 
distributions [2,3]. 

The uncertainties analysis was performed for the CJK photon parametrization only. The set 
of the best values of parameters n,a,/3 and N v , corresponding to the minimal x 2 °f the global 
fit, Xo (table HJ, is denoted as the S° parametrization. Using the Hessian method we created an 
additional basis of the test parametrizations of the CJK parton densities, {S k , k = 1, • • • ,4}, 
where 4 corresponds to the number of free parameters of the model. The set {S^} allows for the 
calculation of the uncertainty of any physical observable X depending on the parton densities. 
Its best value is given as X(S°). The uncertainty of X, for a displacement from the parton 
densities minimum by A% 2 = T 2 (T - the tolerance parameter) can be calculated with a very 
simple expression (named as master equation by the CTEQ Collaboration) 

where parameter t = 5 in that particular case. After a detailed test of the allowed deviation of 
the global fit from the minimum, we found that T should lie in the range 5 ~ 10 [3]. Note that 
having calculated AX for one value of the tolerance parameter T we can obtain the uncertainty 
of X for any other T by simple scaling of AX. This way sets of {S^} parton densities give 
us a perfect tool for studying of the uncertainties of other physical quantities. One of such 
quantities can be the parton densities themselves. 

In Figure El the up-, strange- and charm-quark and gluon densities calculated in the 
FFNScjif models and the GRV LO, GRS LO and SaSID parametrizations are compared with 
the CJK predictions. We plot for Q 2 = 10 and 100 GeV 2 the ratios q 1 (Other model) /q 1 (CJK) 
and g 7 (0ther parametrization) /q 1 { CJK) of the parton q 1 densities calculated in the CJK model 
and other models (or parametrizations). Solid lines show the CJK fit uncertainties for A% 2 = 25 
computed with the {S^} test parametrizations. First we notice that predictions of FFNScjk 
models in the case of all parton distributions lie between the lines of the CJK uncertain- 
ties. There is only one range of x, namely 0.01^x^0.1 at Q 2 = 10 GeV 2 where the up- and 
down-quark densities predicted by the FFNScjjf 1 fit go slightly beyond the uncertainty band. 
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This indicates that the choice of A% 2 = 25 agrees with the differences among our four models. 
Moreover, the GRV LO parametrization predictions are nearly contained within the CJK model 
uncertainties. Obviously, that is not the case, for the heavy-quark densities. The SaSID results 
differ very substantially from the CJK ones. Next, we observe that, as expected, the up-quark 
distribution is the one best constrained by the experimental data, while the greatest uncertain- 
ties are connected with the gluon densities. In the case of w 7 , the A% 2 = 25 band widens in 
the small x region. Alike in the case of other quark uncertainties it shrinks at high x, this is 
due to the large w-quark density in the photon at high x. On contrary the gluon distributions 
are least constrained at the region of x — > 1. Finally we observe that all uncertainties become 
slightly smaller when we go to higher Q 2 from 10 to 100 GeV 2 (not shown). 

5 Summary 

We enlarged and improved our previous analysis [1] . Here we performed 3 new global fits to the 
F^x^Q 2 ) data, excluding the TPC27 experiment. Two additional models were analysed. New 
fits gave x 2 P er degree of freedom, 1.5-1.7, about 0.25 better than the old results. All features 
of the CJKL model, such as the shape of the heavy-quark distributions, good description of 
the LEP data on the Q 2 dependence of the F^ and on F^ are preserved. We checked that the 
gluon densities of our models agree with the HI measurement of the G 1 distribution performed 
at Q 2 = 74 GeV 2 [35]. 

An analysis of the uncertainties of the CJK parton distributions due to the experimental 
errors based on the Hessian method was performed for the very first time for the photon 
[2,3]. We constructed set of test parametrizations for the CJK model. It allows to computate 
uncertainties of any physical quantity depending on the real photon parton densities. 

Fortran parametrization programs for all models, obtained through parametrization of the 
fit results on the grid, can be obtained from the web-page [27]. Set of the data used in the fits 
is also given there. 
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Figure 1: Predictions for the F^(x, Q 2 )/a for the CJK and FFNScj^ models compared with 
the experimental data [14-24], for small and medium Q 2 as a function of x (logarithmic scale). 
If a few values of Q 2 are displayed in the panel, the average of the smallest and biggest Q 2 was 
taken in the computation. 
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Figure 2: The same as in Fig. d for Q 2 > 20GeV 2 . 
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Figure 3: The same as in Fig. ^for a linear scale in x. 
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Figure 4: The same as in Fig. |2]ibr a linear scale in x. 
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Figure 5: Comparison of the recent OPAL data [24] for the Q 2 -dependence of the averaged over 
0.1 < x < 0.6 F2 /a with the predictions of the CJK and FFNScjif models. 
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Figure 6: As in Fig. El for F^x, Q 2 )/a, averaged over four different x ranges. 
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Figure 7: Comparison of the parton densities predicted by various models at Q 2 = 10 GeV 2 , 
as a function of x. 
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Figure 8: Comparison of the up-quark density at four values of Q 2 in the CJK and FFNScjif 
models, as a function of x. 
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Figure 9: The same as in Fig. |HI for the gluon density. 
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Figure 10: The same as in Fig. |HI for the charm-quark density. 
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Figure 11: Comparison of various contributions to the photon structure function F^^x, Q 2 )/a 
in the CJK model for Q 2 = 5, 20, 100 and 1000 GeV 2 . 
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Figure 12: Comparison of the structure function i^ c (x, Q 2 )/a calculated in the CJK and 
FFNScxft- models and in GRS LO [30] and SaSID [33] parametrizations with the OPAL mea- 
surement [34]. 




Figure 13: Comparison of the gluon distribution obtained in the HI measurement performed at 
Q 2 = 74 GeV 2 [35] with the predictions of the CJK, FFNS C j^l & 2 models and GRV LO [5], 
GRS LO [30] and SaSID [33] parametrizations with the OPAL measurement [34]. 
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Figure 14: Parton densities calculated in FFNS C jk models and GRV LO [5], GRS LO [30] and 
SaSID [33] parametrizations compared with the CJK predictions. We plot for Q 2 = 10 GeV 2 
the q 1 (Other parametrization) /g 7 (CJK) ratios of the parton density calculated in the CJK 
model and its values obtained with other models and parametrizations. Solid lines show the 
CJK fit uncertainties for A% 2 = 25 computed with the set of the {S^} test parametrizations. 
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